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ABSTRACT 



An investigation is made of the drift motion of small particles 
under the influence of acoustic oscillations. The investigation is made 
to determine if the motion has a magnitude great enough to produce 
significant changes on the fuel distribution in the chamber of a 
liquid propellant rocket motor. The calculations are made for the 
motion in both a rectangular and a cylindrical chamber. 

In the rectangular chamber the gas oscillation is restricted 
to the fundamental transverse mode and motion in only one dimension 
is considered. The particle drift velocity, that is the non -oscillating, 
non-damped term in the expression for the particle velocity is found 
in the solution of the second order equation. 

For particle motion in the cylindrical enclosure, only gas 
oscillations in the first transverse or sloshing mode is considered 
and motion is restricted to a transverse plane of the cylinder. The 
particle drift velocity, again a second order term, is determined. 

The magnitude of the drift velocity is calculated using condi- 
tions found in a liquid propellant rocket combustion chamber. Dis- 
tances a typical fuel droplet would move during its average life time 
are calculated. The distances are small compared to tlie size of 
most rocket combustion chambers. 
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I. INTRODUCTION 

It is a well known fact that liquid propellant rocket combus- 
tion chambers often experience a phenomenon called unstable com- 
bustion. The gas properties, such as pressure, density and 
temperature, fluctuate during periods of unstable combustion. 

The fluctuations are not random and definite correlation exists 
between the fluctuations at different positions in space or different 
instants in time. Characteristiceilly, the fluctuations have a fixed 
frequency and the amplitude of the variations can grow to values 
comparable to the steady state value of the property. 

The regularity and the large amplitude of the fluctuations 
distinguish the property variations observed during unstable 
combustion from the random variation in gas properties observed 
during "steady" operation of a combustion chamber. 

The observed frequencies of oscillation vary from 10 to 
as high as 10,000 cycles per second and depend on the driving 
mechanism or phenomenon and the chamber size. The frequen- 
cies can be divided into two general categdries: (1) a low fre- 
quency or chugging instability with frequencies in the neighborhood 
of 100 cycles per second and (2) a high frequency or screeching 
instability with frequencies in the range of 1000 to 10,000 cycles 
per second. 

The low frequency phenomena are often caused by an 
interaction between the rates of propellant injection, combustion 
and rate of ejection of material from the chamber. 
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In the high frequency case the variations of gas properties, 
during unstable combustion, are accompanied by motion of the gas 
in the chamber. In general, the motion of tlie gas and variations 
of properties are similar to the conditions which would be gener- 
ated by a strong acoustic oscillation. This observation suggests 
that the large amplitude variations observed during high frequency 
unstable combustion are a result of a resonance which occurs 
when a driving mechanism, associated with the combustion phe- 
nomenon, is in phase with the natural acoustic modes of the com- 
bustion chamber. For example, in cylindrical chambers, the 
oscillations occur in the radial, tangential and longitudinal modes 
plus their various combinations. The first transverse or sloshing 
mode is the mode of oscillation most frequently found experimen- 
tally^^^’^^^ in combustion chambers with high frequency instability. 

The present investigation will be concerned with the high 
frequency instability. A complete and satisfactory answer to the 
question as to why this instability occurs has, as yet, not been 
found. In the attempt to answer the above question, many mech- 
anisms have been proposed that are based on the interaction of the 
pressure atnd temperature oscillations with the combustion process. 
The interaction depends on the existence of a time delay which is 
essentially the delay time in the response of the burning rate 
to a change in temperature or pressure. 

Since the delay times in combustion systems are very 
sensitive functions of the fuel-oxidizer mixture ratios, a systematic 
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variation of the mb:turc ratio in the chaxnber may set up conditions 
tliat are favorable for unstable burning. The suggestion has been 
made that variations in the fuel>oxidizer ratio can be produced by 
the relative mean motion of droplets of fuel and oxidizer under 
the influence of acoustic oscillations in the combustion chamber 
gas. 

This investigation is concerned with the initial require* 
ments for the production of the mean particle motion mentioned 
above. The requirements are first, that a particle drift motion 
does occur in a chamber with acoustic oscillations and, second, 
that the magnitude be large enough to result in the appreciable 
particle motion which would be necessary to produce significant 
changes in fuel composition in the chamber. 

The following two idealized problems will be solved. 

First, as a simple example, particle motion in a rectangular 
chamber with acoustic oscillations in tlie fundamental transverse 
mode is considered. In the second, more complex problem the 
calculation is carried out assuming a cylindrical enclosure with 
acoustic oscillations. The small particle drift velocity induced 
by the first transverse mode of oscillation is determined. 

To calculate the particle drift velocity the equation for 
the velocity of the gas undergoing acoustic oscillations in the 
chamber is combined with Stokes drag law to calculate the forces 
acting on small particles. The resulting equation is linearized 
by introducing a perturbation expansion in powers of a term 
proportional to the amplitude of the drag force acting on the 
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particle as a result of the gas motion. It io necessary to go to 
second order terms in the solution for the particle position in 
order to obtain a mean drift velocity. 

Representative values of the drift velocities are calculated 
using conditions found in a liquid propellant rocket combustion 
chamber. Directions and distsuices the droplets would move, in 
their average life time, are determined. Finally, an upper limit 
is obtained for a maximum variation in particle density as a re- 
sult of the drift velocity. 
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II. RECTANGULAR CHAMBER 



Development of the Equations. Before proceeding to the more 
complex and more realistic problem of the investigation of par- 
ticle motion in a cylindrical enclosure with acoustic oscillations 
it will be illuminating to look at particle motion in a rectangular 
enclosure containing acoustic waves. The rectangular case will 
illustrate the important features of the problem in a straight- 
forward manner since motion in just one dimension need be con- 
sidered. 

The method of solution for the particle drift velocity will 
be as outlined in the introduction. 

In solving the problem for the case of the rectangular 
chamber the following idealized assumptions will be made: 

1. The fluid medium is considered to be at rest and in a uni- 
form state e>:cept for weak isoenergetic fluctuations in velocity 
and state propearties. 

2. The small particles are considered to be spherical and to 
experience Stokes law forces. 

3. The particles have negligible effect upon the gas motion in 
the chamber. 

Consider a rectangular chamber, as shown below, with 
acoustic oscillations in the fundamental transverse mode. 
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Gas motion io restricted to oscillations in the direction. The 
linearized partial differential equation for the gas pressure dis- 
turbance in the chamber, under the conditions considered here, 
is a simple wave equation of the form: 



^ ^ _L ^ 



( 2 . 1 ) 



The well-known solution to equation 2. 1, restricting the solution 
to real values, is: 



P = P sin CoS cot 

' I o n 



( 2 . 2 ) 



Here, R is the amplitude of pressure fluctuations of the gas in 
the chamber, C is the local velocity of sound, CO is ZTf times 
the frequency, , and t is the time. CO must be to sat- 

isfy the boundary condition that the gas velocity and, hence, the 
pressure gradient vanish at the chamber v/alls ('y. =-^ ). For the 
first transverse mode considered here ^ is 1. The velocity and 
pressure fluctuations of the gas are related by the following equa- 
tion: 



:= - VP = 

The velocity of the gas in the chamber is given by integration 
of 2.3 as: 

U = Cos ^ sin cot , 



(2.3) 



(2.4) 
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Equation 2. 4 shovrs that the gas velocity is zero at the chamber 
walls and a maKimum at the center of the chamber. The ampli- 
tude of the gas velocity oscillations is proportional to cosJ5^ . 

If the location of the particle is given by f , and if assump- 
tion 2 holds, then the drag force on the particle is given by Stokes 
law and is: 



The equation of motion for the particle is obtained by equating tlie 
inertial forces and drag force, and is: 



It is important to note that the gas velocity must be evaluated at 
the position of the particle, that iQ U = U(fjt) • M assumption 
3 holds, the gas speed, U , can be eliminated by substituting equa- 
tion 2.4 into 2.6; this gives: 



This result can be simplified by the use of the following dimen- 
sionless variables: 



Introducing the dimensionless variables in equation 2.7 gives: 




(2.5) 






( 2 . 8 ) 



-f- ~ t^'cos X ^ sin sirD" 

jy 



(2.9) 




% 
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Here, the constants are defined as: 




rm c 



9nd 




- ^6 TT/Ze-' j? H 



( 2 . 10 ) 



The solution to equation 2. 9 will give the position of the small 
particles, ^ , at any time.O^ , while they are under the influ- 
ence of the oscillating gas in the chamber. 

Note that equation 2.9 is non-linear due to the appearance 
of 2 in the trigonometric function on the right hand side of the 
equation. Because of this nonlinearity, an approximate method 

of solution of the equation must be used. The approximation used 

/ 

depends on the fact that the term £ can be held to small values 

by limiting the magnitude of the viscosity, or other flow proper- 

/ 

ties appearing in S . Note that this restriction on the calculation 
is Independent of that required to obtain a linear acoustic equa- 
tion (2. 1) in the first place. Thus, when the expansion in terms 
of £ is continued to terms of higher order than the first, there 
is no obligation to extend the approximation in the acoustic equa- 
tion to higher powers of pressure or velocity amplitude. By re- 

/ 

stricting the analysis to arbitrarily small values of £ , it is 
possible to linearize equation 2.9 by introducing a perturbation 
expansion of the form: 

?(« = zjo-) (-y) + . {2.U) 

Substituting the expansion into equation 2. 9 and separating the 
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equation into equations for the various orders of £ 



gives: 







0 



( 2 . 12 ) 



-f 1/ -d^ - CoS ^ io sin (2.13) 



S/'na-TrT . 

In making the above substitution, Cos is expanded as 



(2. 14) 






For a particle at rest at 0^= O and at the point 2 » 

appropriate initial conditions are: 

= 2 ^ = 0 
?, (o') = ^(0) = O 

( 0 ) = ^ ^ ■ 

Solutions to equations 2. 12, 2. 13, and 2. 14 contain terms 
-which are constant, oscillatory, exponential and linear functions 
of the time. The purely oscillatory solutions to the above equa- 
tions are not of interest since no net transport of particles is 
involved and the solutions involving decaying exponentials have 
a negligible effect on the drift velocity after a few oscillations. 
Hence, the only parts of the solution that are of interest and that 
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result in a particle drift velocity are those solutions which are 
proportional to time . 

LfOoking at equatioiis 2. 12 and 2.13, it is seen that the so- 
lutions to these equations are of the form of constants, decaying 
exponentials, and sinusoidal terms. However, equation 2. 14 
has Z, , multiplying the right hand side so that the product of the 
trigonometric solutions of Z, with [sinS'frD'J will prochace a 
constant term. The solution for will then contain the first 
and presumably most important term proportional to time. 

in the following paragraphs the solutions for , Z, 
and ^ will be found. The chief aim of the analysis is to ob- 
tain the single term which is proportional to time so that tlie 
drift velocity, as defined here, can be evaluated. 

Solutions to the Squations. By inspection it is obvious that the 
solution to equation 2.12, with the given boundary conditions, 
ie ^ - constant. 



Substituting Zo = Z in equation 2. 13 leads to: 




( 2 . 16 ) 



The solution to the homogeneous equation is: 



2, = C, * C^e 



(2.17) 



The trial particular solution is: 



Z = C, sin 3.rry + C. cos arrCT 

I 3 ^ 



(2.18) 
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Substituting equation 2. 18 into 2. 16 gives: 

- sin s.'rZr - 4'7r^Qcos -f- r cos 5rT 

^l^'sTrC^ sin S.^ZT' - jcos^^ S/'n^Tr'y 



(2.19) 



Equating the coefficients of like terms gives the following values 
for the constants: 



r - 

^3 -_4T^ 



Q - K ^ 
^ -KlsTr-'^ir^ 



( 2 . 20 ) 



The solution for ?, can now be written as: 

-7 ^ ^ CoS^^ -V k 'cos^^ ^ 

Z, - C, -/-u e -f — 5-*^ SmsTj-h ^—3 

L -4ir=*-K^ -k^-Tr-^ir^ 

The initial conditions that 2, ("(J) = ~^(o) = Q are used to evzd- 

a'U 

uate the constants C, eind . The following expression is 
obtained for the complete solution to equation 2. 13: 



( 2 . 21 ) 



2 = 



3TC0S^^ 1<C0S%^ 



^Tt^K'+ K'' 



, _ kiy 

-‘^'rkk'- k'^ ^ 



-j- 5/V?5^Cr k QOS SlTT''^ 

-4Tr^--is'^ -^^*57r'-S? 7r^ 



■K^^TT-^TT 



(2.22) 



The differential equation for 2^ can be found by sub- 
stituting 2. 22 and =2 into equation 2. 14: 



c/ -^a , r 









^ \-4ir^K'- k‘ 



^L.' u'3 J <S 



S/n. siTT 



r 



- i sin it (s/o 

Sia^Try cos ^rr'y 

^ V-/c^r-^7rV 



(2.23) 






..Jii 
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The above equation may be simplified slightly by using the follov/- 
ing trigonometric identities: 



sin cos - 

Sin srT cos JJttT = 
(s/n5 7rCr/= 



^ Sin irZ. 



sin. ■4- ?r iy 



cos ^TT CT 



(2.24) 



Equation 2. 23 reduces to: 



^ ciy 



sin 



^ r-'' - ^ 

< 27 rj -he sin^yrO^ 




(2.25) 



Here, 



D = 

z' 

F = 



'TTsin. rZ. 
n^si'n ir^ 



n~ — ^ S m Tt " k 

^ " k {^rr^k'^) '■ 

C ' = /r s in '7C:i_ 



-^Ar'(4/kzk'^J ' 

The solution to the homogeneous equation is again: 

cz r r - 
^ - C, y- Q e 

The following trial solution will be used to find the particular solu- 
tion to equation 2. 25. 

^2 ~ [fj'^ sin 5 ^ if -h COS SLTr CT 

^ -K'a . ^ ^ -Fy 

-i-Q C 5//2 5 '^J -z- C-p e. cos a ttJ 



Cy Qos^-rr'Cy -t Ccj sin. J 



(2.26) 
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Substituting 2. 26 Into 2. 25, equating coefficients of like terms, 
and solving for the constants gives: 



- 



C, + Q e 






nr sin, ir'2 



r-f, 






jcrs^rX 



-K'd' 



+ 



^ — ■ / 



^ ^ , ]sfn4nry -r 



/'^7r(Fi-&') 



sin 3.7T J"' 



(2.27) 



rk'F' _ 'i . /rtr 

Here it is unnecessary to evaluate the constants C, and 
since C, is a constant and will disappear when the derivative 
v/ith respect to time is taken to find the particle velocity and 
is multiplying a decaying exponential which will not contribute 
to the steady state particle drift velocity. The non-oscillating, 
non-decaying term, a second order term, is ^ ^'')J 

Using equation 2.11, the non-oscillating and non-decaying part 
cf the expression for the particle position can be written as: 




z - 



'EF/n TT^ _ 



yP 



(2.28) 



Or, using 2.3, equation 2. 28 can be written in terms of the origi- 
nal variables as: 






(2.29) 






m 
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Now taking the derivative of 2. 29 with respect to time , the mean 
particle drift velocity is obtained: 






si. 












(2.30) 



Remember here that E and |( = • 

Equation 2. 30 gives the steady state drift velocity of the 
particle which is at a point ^ in the chamber. The drift term 



is second order in E and depends sinusoidally on position. The 
greatest drift velocity occurs at an initial particle position of 




The fact that the drift velocity is second order can be jus- 
tified by the following argument. The drift velocity exists because 
of the gradient in the amplitude of the gas oscillations. Consider 
the forces acting on a particle due to the difference between par- 
ticle and gas velocity. The gas is oscillating at a fixed frequency 
but with an amplitude that varies spatially. Thus, during a 
single oscillation, the particle moves from a region in which the 
amplitude of gas motion is large to a region in which the amplitude 
of gas motion is sJ^Hler. During that part of an oscillation when 
the particle is near its maximum excursion in the direction of in- 
creasing amplitude of gas motion forces act to accelerate the 
particle toward the direction of decreasing amplitude. Similarly, 
when the particle is on the opposite side of its mo 2 ui position, 
forces act to accelerate the particle in the direction of increasing 
amplitude of gas motion. However, since the force acting on the 
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particle is proportional to the difference between the gas speed 
and the particle speed, tlie forces acting on the particle are higher 
in the regions of higher amplitude of gas motion. Konce, over a 
number of oscillations of the particle, a net force acts to accel- 
erate the particles toward the direction of decreasing amplitude 
of gas oscillations. The net driving force depends on the differ- 
ence between the forces acting on the paarticle on either side of 
the mean particle position and thus depends on the mean difference 
in the gas speed on either side of the particle mean position. 
Hence, the driving force and the resulting drift velocity is pro- 
portional to the product of the gradient in amplitude of gas motion 
and the amplitude of particle motion. Since both of these quan- 
tities are first order, the resulting particle drift speed is second 
order. 

The gradient in the amplitude of gas oscillation is pro- 
portional to sin^Z &nd the amplitude of the particle oscillation 
is proportional to cos-^? , e.g. see equation 2.22. Thus, the 
maximum drift speed would be expected at the value of Z which 
maximizes (cos^H)(sin-^z) or ^ * 

Discussion of Solutions. It Is of interest to determine the man- 
ner in which this drift depends upon the p>article size. To do 



this, it is convenient to use the fact that S can be written as a 



nated from equation 2. 30 to give the following expression for the 



/ 



function of 




Therefore, £ ^ can be elimi- 



drift velocity: 
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nr = 







(2.31) 



To make the dependence upon particle radius explicit, write tlie 
particle mass as * In these 

terms the drift velocity is given by: 



nr = 





(2.32) 



As the particle radius , (T' , gets large it is evident from equation 
2. 32 that the drift velocity approaches zero. If equation 2. 32 is 
rewritten as: 



(V - 




(2.33) 



it is evident that as the particle radius gets small, the drift ve- 
locity again approaches zero. Therefore, a maximum drift 
velocity exists for some finite value of particle radius. This op- 
timum radius. o~ , resulting in maximum drift velocity, can be 
found by examining the first and second derivatives with respect to 
O'. The result of the calculation gives: 



■K- 



c r 



(2.34) 



Substituting the optimum particle radius, 2. 34, into the 



drift velocity expression, 2. 32, gives the following expression 
for the maximum drift velocity; 




I 






In terms of the chamber pressure, fc , the ma^imuin drift is: 

Note that while the optimum particle radius is a function of 
the chamber dimension, , the maximum drift velocity, ri/\ is in- 
dependent of the chamber size. Therefore, the drift motion is less 
important as the chamber size increases. 

Note also that the maximum drift velocity is independent of 
droplet properties and only depends on the gas properties through 
the speed of sound C , specific heat ratio Ja and pressure ratio 

t ‘ 

The ratio of drift velocity for arbitrary particle radius to 
the drift velocity for optimum particle radius is simply obtained 
from equations 2. 32 and 2. 34 and is given by 



^ = a ] 

L I 



(2.37) 



The drift speed drops rapidly for particle radius different from 
optimum radius. For example, for 6^ ~ ^ <5" or 3.6~ the ratio 
is about 4 . 

(V * 
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m. CYLINDRICAL CHAMBER 

Development of the Equations. The techniques used in Section II 
are now applied to the solution of the problem of particle motion 
in a cylindrical geometry. The problems are quite similar except 
that with the cylindrical geometry, (A) motion in two directions 
must be considered and (B) centrifugal forces may be important. 

To investigate the motion of small particles in a cylindri- 
cal enclosure under the influence of acoustic oscillations, the 
same assumptions that were made for the rectangular enclosure 
will be used; namely: 

1. The fluid medium is considered to be at rest and in a uniform 
state except for v/eak isoenergetic fluctuations in velocity and 
state properties. 

2. The small particles are spherical and experience Stokes law 
forces. 

3. The particles have negligible effect upon the gas motion in 
the chamber. 

For the conditions given in assvunption 1, it can bo shown 
(3) that the gas motion is governed by the wave equation which is 
v/ritten below in cylindrical coordinates for the pressure distur- 
bance: 



-L IB- 







4 - 



J_ ^ _ _L & 



( 3 . 1 ) 



Since we are interested in finding a particle drift velocity 
in a cross sectional plane of the cylinder, in analogy to the simple 
transverse waves considered in the rectangular chamber, only 
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oscillations independent of 2 will be considered in this analysis. 
The solution to equation 3. 1 for a circular cylinder with stamding 
acoustic oscillations is given by: 

P = (3.; 



Here, as in the previous problem, the frequency, P , is , 

and, in this case the angular velocity, 6^;^, , is 

arise from the condition of zero radial velocity at the cyl- 
inder walls. The subscripts rm and ^ specify the modes of vibra- 
tion in the and /I directions. At the walls is a solution 

to 



— Km (, 



For each value oirrn there exists an infinite number of solutions 
to equation 3. 3 since IXyn Is an oscillatory function of its argument. 
Characteristic vetlues of are listed in reference 3. 

For oscillations restricted to a transverse plane of the 
cylinder the simplest cases are the first radial mode C , 

1. 22), and the first transverse or sloshing mode {rrn- i ,rri~c , 
*586). In investigating the possible particle drift velocities 
induced by the acoustic oscillations only the effect of the sloshing 
mode will be considered since it is this mode which is most often 
observed in liquid and solid propellant rocket systems experiencing 
high frequency instability. 

With the restrictions to oscillations independent of 2 ! and 
to real solutions, equation 3.2 can be written: 



P = g cos f 3", cos 37TVT . 



<3.4) 



The coordinate system used in this problem is shown in the follow- 
ing diagram: 




In this case: 

+ i [f?(«3c^)(c«5^.t)(5J.”)[ 74 "^")- ( 3 -^'^)]} _ < 3 - 5 ) 

An expression for the gas velocity can be obtained from equation 
3. 4 in a mamner analogous to that used in the case of the rectangu- 
lar chamber. The result of the calculation is: 

0 = [sM^ sm if y, s/n ^Trpt] 



-f i 



Po 



Vi./ 4\)^a 



cos^ sMsr^T (3.6) 



Fig. i shows the magnitude and direction of the gas velocity, 
as a function of , for various points in the cross section of 

the cylinder. The calculations are made for a time such that 
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sin ^'Tr\>t = 1. Since the cross section is symmetrical, only the 

first two quadrants are shown in the plot. It can be seen in Fig. 

o o 

1 that there is a gas velocity node at (f equals 0 and 180 while 
the pressure node occurs at (p = 90® and 270®. 

Using assumption 2 the drag force on the particles is given 

by: 



^ 7 T^ ^ (^0 



(3.7) 



Here, as before, U is the gas speed and nr is the particle speed. 
The particle position is given by /i , the particle velocity by 

, and the particle acceleration by 



-4- h 



fJi - t Ua -a/m 
7t - (at 



+ (-tn 



a 



M- -/. /I d^(f’ 1 



JT d't 



dt 






The equation of motion for the particle can be obtained by equating 
the particle inertial forces to the drag force obtained from Stokes 
law. The result is: 








rrn ^ 



d'l djL 

'-at dt 







(3.8) 



-tL 






In the above equation the same symbols are used for gas position 
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as for particle position since the gas velocity must be evaluated 
at the position of the particle (by assumption 3). Equation 3.8 
can be separated into two more useful equations corresponding to 
the and directions: 

Ln, component 



6 \ , (>iryU<s- (J/i. _ — 

^ [dtj - 

fo Cos siri 

5^ ^ Vq. I L o\ a- / -av a JJ 

Lt^ component 



(3.9) 





(3.40) 



As in the case of the rectangular chamber it will be conven- 
ient to nondimensionalisse equations 3.9 and 3.10; the following 
dimensionless variables are used: 

i = ^ ‘f ■ (3- “) 

Introducing the dimensionless variables into equation 3.9 and 3. 10 
gives: 

component 






+ - 




(3.12) 



- £ CoS (|) sin 



(7rcx;,2) - ir^(7ro^,o 
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■ (f component 



jrjr ^ ^ Jr^ 

S sin (p Sin 3; J 



(3.13) 



Here, g= 3//er-^ 2 m^ are parameters of the prob- 

^ /7770 yc* /^c ^ 

/ / 

lem and, except for constants, are the same as the & and K ha 
the rectangular chamber problem. 

Equations 3. 12 and 3.13 are non-linear due to the appear- 
ance of ? in the argument of the Bessel functions cuid are analo- 
gous to equation 2.9 in the case of the rectangular chamber. Since 
E is proportional to fl and other flov/ properties, an argu- 

ment similar to that made about equation 2. 9 is valid. The follow- 
ing perturbation expansion in powers of E is introduced to linearise 
equations 3. 12 auid 3. 13. 



= HjTj + e +--- <3.i4) 

(+>(•7) = -fSi, (r) ^ + — o.is) 



In substituting the above expansions in equations 3. 12 and 3. 13, 
sin (p is expanded as S/nf = (s/n p^)(i-€^f(p^^e(p,r-h ---) 

+ (coS<pt,)(€ 4’i } and cos is expanded as 

COS(p = - (s/n(pj(t(f, i-S%- . 

The Bessel functions are expzinded in a Taylor series about the point 
to get: 
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'X -ez, X 2 o) 7r<x,, 






(-r=<,f T. (Tr«,.Z.) + 3;f-7r«’„Zj] 



T,('^^,o^) = fx ?o) -h s T, (T^roX ) 

^ ITo 

+£V=(„?^3;(ip=<.z„) - -^ 3; fir«„z,) -^-^%-<,.f[p^f~']y,(Tr^„iO 



U^ ( TT <^,0 Z) = 



5 (’>'-io ?J ■^ f Z, ir<^,. 3 : ( •»•«„ ?, ) - £H' 3Q (t«„Z.) 
-f- X (rc^„?o) - 

*^o 5. 

3;(7r=s.?„) + 



The results of substituting 3. 14 and 3. 15 into equations 
3.12 and 3. 13 and separating the results into orders of £ are: 

L/l component 

4 . ^ (dis f' ~ 0 (3.16) 



IIl 

dcr 



u di> -3^ d^i 

^ dr "dr dr 



' (#/ = 



fe- s>'^ cos (f>. -t- j; I 

1 ^ 



(3.17) 
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p- 



sm 7r<=<,^ cos (f)^ ZT, [Tf^,o ^o) 



- ^ cos + <f>, sin 4t,^X('^-^,.^,) -XCir‘^,.s.jl < 5 -‘ 8 ) 



/ (P component 



sz, 4^ ^ -f -i- ?/ k 4% - 0 



jr dr d j 



ar 



(3.19) 



^^“drdr ^drdcr 'dr c/r dr^ 



dr* 



+ 2 .'^^ + 5 K?.^ ^ = 



Sinp^ sin a^9J j (3.20) 



d\ 



-^‘'^4 fSZ^i,d% +p^ 



dr’ 



dr“ 



^ ' ^r- 



4^4^! +-<^r2 

O' r 



^r 



4 v ^ ^ - i^' -^o 1'^^'° ^ - i- - . r 

c^cr c>2T t ‘- ^ ' 

-h (^, cos p^ J, sin . 

If a particle is considered to be at rest at C - C 



0l (3. 21) 

and at a 



point 2^ y 4' . then the appropriate initial conditions are: 



Zo t^) - H j <pj^) ^ 4>; g^{»)^ ( 0 ) ^ c '] 



H, (0) = 


d ^1 ( ^ ^ _ 
’ dr '' 


<p, 


(f'} 


- ^ l\ (0) =: C 

CJ 


s 

( 

\ 


(3.22) 




zr ^ (S) = 
S3 ^ 


>ti' 


Xj 


=: vO^ = C 

J 


) 
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It i3 noted again that purely constant, damped or oscillatory 
solutions will not be o£ interest since no net transport of particles 
is involved. 

The form of equations 3. 16 to 3. 21 cam be compared with 
equations 2. 12 to 2.14. Since (j)^ , and multiply the right side 
of the second order equations it is expected that any terms in the 
solution which are proportional to time will again appear in the sec- 
ond order equations 3. 18 and 3. 21. 

The expectation that the drift velocity will be second order 
in £ , specific £illy for the radial motion of the particles, also fol- 
lows from the physical arguments given on page 14. However, 
for the geometry considered here, centrifugal forces also act on 
the particle. Since the centrifugal forces are proportional to the 
square of the particle velocity and the particle velocity is first 
order in £ , the drift which results from centrifugal forces is 
also a second older term* 

Solutions to the Equations. A solution to the sero order equations 
3. 16 and 3. 19 is = constamt = 2 and <P^ = constant = (p . The 
constant £ is closely related to a drag force and can be called 
2 m interaction parameter so that to zero order in £ there is no 
drag or interaction between the particles and the gas in the cham- 
ber with the result that the particle position is fixed. That is to 
zero order in S , the particles are not moved by the gas, but are 
at a fixed position. Using these simple zero order solutions, the 
solutions for , 2^ * 0, , and can be found by the same 
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method used in solving equations Z. 13 and 2. 14 la the rectangular 
chaniber problem. The complete solutions to equations 3.17, 3. 18, 
3. ZO and 3. Z1 are discussed in more detail in Appendix A. 

The first and presumably most important non-oscillatory 
and non-damped terms appear in the solutions to the second order 
equations 3. 18 and 3. 21. These non-oscillatory and non-decaying 
solutions are: 



Here, 



D 

E 

F 

G 

H 

X 



^(F£ + D6)l ^ 

^ f irSl" tSil! r 

r 1^/.? ? [t^ ^) - 2 ) j 

r cos (J [(7r-<, , ( rcx',^Z) - U' ( )] 

= sm ^[3; - T^(ir<=-,J)] 

- Sin $ [r^,o I 0 ('TT«/o2) ^ 

- Cos ^ HT 



(3.23) 



Differentiation of equation 3. 23 with respect to time will 
give the particle drift velocity. The objection might be made that 
the decaying exponential terms in the solutions to the first order 
equations may contribute appreciably to the particle motion. A 






i i 
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calculation to determine the magnitude o£ the first order decaying 
exponential terms, compared with the steady state drift velocity, 
is shown in Appendix B. In general, these terms appear to be neg- 
ligible for reasonable values of the parameters. 

Moreover, the terms which decay exponentially in time 
arise from certain assumed initial conditions of the particle posi- 
tion in a certain phase of the motion. Since there is no physical 
way in which one initial condition or phase may be preferred over 
another, these can occur at random and hence the decaying terms 
can not contribute to the mean motion of large masses of particles. 



IV. RESULTS AND DISCUSSION 



The Drift Velocity. Ignoring the oscillatory terms and decaying 
exponentials the drift velocity of the particles can now be written 
down using the solutions listed in 3. 23. Remember that the velocity 
of the particles is 

or in the dimensionless expanded form is: 



+ £ ^ t +£“ 2 , 6=2, A 



The component of the drift velocity can be written as: 



(4.1) 






- ^ r ? D^~ -tDQ) ! £^c r ~^ - (F E+ D&) I ,A 7\ 



in terms of the original variables. Note that (Ft +D&)] 

is a constant for atny given /[ , (p and cylinder radius and does 
not depend on the particle dimensions. Let (ft'-hDG)] - ^ j 

then the JZ component of drift velocity is more conveniently 
written as 



_ e“c<? 



(4. 3) 



To see how the A, component of particle velocity varies with par- 
ticle size the expressions for C cUid )<( can be substituted into 
; this substitution gives: 
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rvx - 



16 



,1^)^ 

[[v^J 



I 



Or, since the mass of the particle equals 
be written as 






J_ Mff3. 

S’/, 

!6> ^pTT^)>^<r'^ 



(4.4) 



4. 3 can 



(4.5) 



As a~'-^oo in equation 4. 5, C> and as cT'— >■ O , ^ Q ; 

therefore, as in the case of the rectangular chamber, there is an 
optimum particle radius, cr , which results in a maximum ^ • 
This optimum radius is again found by examining the first and 
second derivatives of the velocity with respect to 6 " . The result 
of the calculsition gives: 



6~' 

A 





(4.6) 



The same arguments can also be made for the component. 
Thus, the (f component of particle velocity is 

nr - -Si f~^ +0^) I ^ f j 

^ I /r/> V - /a'/* i 



(4.7) 



in terms of the original variables. Again, ^ (Bh+Cl) is constant 
for a given Ti , (f and cylinder radius. If -^u { £h^0t) - ^ ^ 
then the ^ component of drift velocity can be written as 







\ c ^ 



(4.8) 






1 / 
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Since the expression for the A, component o£ drift velocity (equa- 
tion 4. 3) differs from the cp component of drift velocity (equation 
4.8) only by a multiplying constant, it is evident that the optimum 

radius for maximum will again be . 

T if? V4p^V7r 

Since the optimum particle radius is the same for the A- 
and (f components of drift velocity, a particle of radius <r will 
indeed experience a maximum resultant drift velocity. Note that, 
since ^ in the case of the rectangular chamber, the optimum 
radius given in 4. 6 agrees exactly with the optimum radius found 
in equation 2.34 for particle motion in a rectangular chamber. 

Again, in analogy with the rectangular chamber problem, 
it will be interesting to see whether or not the maximum particle 
drift velocity is independent of the cylinder radius, CL . The par- 
ticle drift velocity can be written as: 



a/' = i. 






_ f t ^ ^ 



(4.9) 



i 



h 



Here ^ and S contain the cylinder radius in the ratio A.i and 

for a given they are independent of the radius. It remains to 

be shown that S , K and are independent of the cylinder 

radius. Since "0= , it is obvious that V’ varies as the In- 

di- 
verse of the radius, and that does not vary with the 

radius. Since the optimum particle radius varies as a* and the 

mass of an optimum sized particle varies as (2 it is evident that 

r = and Lc= 6 u are both Independent of the cylinder 

radius. Therefore, while the optimum particle radius is a function 



of the cylinder size, the magnitude of the maximum particle drift 
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velocity is not. This result is in agreement with the result ob- 
tained for the rectangular chamber problem. 

Another comparison with the rectangular chamber problem 
can now be made. If the optimum particle radius, equation 4.6, is 
substituted into the particle drift velocity expression, equation 4.9, 
the following is obtained: 



- 1 . 

nr - 



c 

4^"^ (7r=x,o)^ 




(4. 10) 



For Cp = 0 > the constant is Q and the particle motion is en- 
tirely in the radial direction. Fig. 1 shows that the gas velocity 

O 

in the cylinder is also entirely in the radial direction when (p is O . 
Since motion was restricted to one dimension ih the rectangular geometry 
problem, a comparison of the maximum particle drift velocities is 

O 



most meaningful if compared for a cp of 0 in the cylindrical problem. 

O 

When cp = O is substituted into equation 4. 10, the equation reduces to: 



For comparison, equation 2.36 is rewritten here: 

\SL 



X 

nr - 



sin. 



(2.36) 



The expressions are the same except for the terms on the right re- 
sulting from the geometry difference in the two problems. 

It is also noted that the ratio of particle drift velocity for 
arbitrary radius to the drift velocity for particles of optimum radius 
varies in exactly the same manner as given in equation 2. 37 for the 
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rectangular chamber problem. 

Magnitude of the Drift Velocity and Discussion. Thus far, with the 
assumptions listed in Section 11 and the restrictions mentioned for 
S, the expressions for particle drift velocity are perfectly general. 
The expression for particle drift velocity, equation 4.9, will now 
be used to calculate the drift velocity of small liquid droplets in a 
liquid propellant rocket combustion chamber. Although the fluid 
medium in the combustion chamber is not at rest, the Mach number 
of the gas relative to the particle is usually low near the injector. 
Stokes law drag will hold if the droplet shape remains nearly spher- 
ical and if the Reynolds number of the flow is less than one which, 
for reasonable values of the parameters, means that the particle 
radius must be less than ten microns. The results of these drift 
velocity calculations will not be exact, but will give an indication 
of the mean motion the droplets experience in a combustion cham- 
ber. 

The follo'vving properties, for the droplets and gas in the 
chamber will be assumed: 

For the particle: 

= 2. 1 slug/ft. ^ 

The particle radius will be picked to 

maximize the drift velocity amd for 

the values listed here: 

-5 

6”' = 2.75 X lO” ft. =8.4 microns 
For the gas in the chamber: 

C = 2360 ft. / sec . 

= 1.54x i0'° alug/sec.ft. 

P = 4. 46 X 10*^ alucr/ft. ^ 
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J. =1.3 
T = 2500® K 
P = 692 1/ sec. 

For the chamber: 



a = 1 ft. 

In substituting the above values in the equations the form of £ 
was chauiged to £ = ^ substituting ^ ^ for C* • 

Here, and ^ are the mean chamber pressure and density. 

In calculating the magnitudes of the drift velocities, values of 
A of . 1, .5, and . 9 ft. were chosen for the initied particle posi- 
tion, and since the drift velocities are completely symmetriced, 
values of (f in the first two quadrants of the cylinder cross sec- 
tion were considered. In particular, angles of 0®, 22.5®, 45®, 
67.5®, 90®, and the corresponding angles in tlie second quadrant 
were used. Table i shows the results of the drift velocity calcu- 
lations. In both Table I and Fig. 2 the magnitude of the particle 
drift velocity is divided by • Sizice the amplitude of 

pressure fluctuations in an unstable rocket chamber is a parameter 
vdiich varies quite drastically, it is convenient to multiply the values 
listed in Table I and Fig. 2 by the square of any desired ratio of 
pressure amplitude to mean chamber pressure to determine the 
actual drift velocity. 

Fig. 2 shov/s the directions and magnitudes of the drift 
velocities in the first two quadrants of the cylinder cross section. 

It is seen that the drift near the wall of the cylinder is greatest at 
= 90® which is to be expected since at this value of (f there 
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is a pressure node and, as cajQ be seen in Fig. 1, the gas velocity 
is entirely in the (f> direction. The drift velocities are generally 
in the radial direction, due to the centrifugal force exerted on the 
particles by the oscillating gas. The drift velocities are greatest 
when the initial particle position is at the mid-radius of the cylinder. 
Since the drift velocity is a second order phenomenon, it is difficult 
to find any correlation between the drift direction and magnitude, 
and the pressure and velocity distribution of the oscillating gas in 
the chamber except for that mentioned above. 



amplitude amounting to 31.6 per cent of tlie mean chamber pres- 



I, the particle drift velocity varies from . 39 to 4 feet per second. 

As mentioned earlier, the magnitude of the first order de- 
caying exponential terms in the solution of the cylindrical chamber 
problem are investigated in Appendix B. For tiiis numerical 
example, the results show that after just two oscillations the mag- 
nitude of the velocity contribution of the decaying exponential 
terms is leas than one-one hundredth of the magnitude of the second 
order drift velocity. 



To obtain a range of values for the magnitude of the drift 
velocity, a value of . 316 is chosen for ; that is, a pressure 




sure. Then 




is 0. 1 and with the values shown in Table 



It is also interesting to determine how far the particles 



move in the combustion chamber. In .01 seconds, which is close 
to an upper limit for liquid droplet life times in a rocket chamber 
(4), the droplets can move from . 05 to .43 inches depending on their 
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initial position in the cylinder. This is a very small distance when 
compared with the size of most rocket chambers. 

In comparison, the drift velocity in the rectangular chamber 
will reach a maximum of 17.4 feet per second at an initial position 
in the center of the chamber. In . 01 seconds the particles could 
move 2 inches which is considerably greater than in the cylindrical 
chamber, but still not significant \vhen compared with chamber size 
or wave length of the oscillation. 

Variation in Particle Density. Although it is not the aim of this 
present work to calculate in detail the variation in fuel-oxidizer 
ratio which results from the relative drifts of fuel and oxidizer par- 
ticles, it is of interest to calculate a few representative values for 
the change of particle density with time. The variation in particle 
density, . cam be computed directly from the continuity equation 
for nonsteady flows. The equation, in vector form, is 

-t = 0 . (4.12) 

^ is the particle drift velocity is the droplet density 

in number per unit volume. If the problem is restricted to the cal- 
culation of variations in density over short periods, then the varia- 
tion of Ot with time may be approximated by 

(p) * «•! . 

Also, for short periods, the effects of evaporation of the particles 
may be neglected and hence the particle radius will be constant. 

V/ith the above restrictions equation 4. 12 reduces to 
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• (4.13) 

A/>V' 

Values of are calculated for the chamber conditions given 

in the previous section. The calculation is carried out for the cyl- 
indrical chamber at ^ =0.1 euad Cp = 0 and for the rectangular 
chamber at ^ = 0. 1. A time of about 1/5 of the particle life time 
of 0.01 seconds is assumed and a radius equal to the optimum value 
for particle motion is used in the calculation. 

A 3 per cent reduction in particle concentration is obtained 
for the cylindrical chamber and a 10 per cent reduction for the 
rectangular chamber. These values of density variation are the 
maximum which can be obtained for tlie chamber conditions used. 
For the cylindrical chamber, the chamges in particle density appear 
to be too small to cause a large change in oxidizer -fuel ratio and 
in fact the variations of the order of 5 per cent are much less than 
the variations which occur because of the injector design. 
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V. CONCLUSIONS 



It has been demonstrated that small particles in an oscil- 
lating gas field experience a steady state drift velocity in both 
rectangular and cylindrical enclosures. This result is based on 
the assiunptions that (1) the fluid medium is at rest and in a uni- 
form state except for weak isoenergetic fluctuations in velocity and 
state properties, (2) the particles are spherical and experience 
Stokes law forces, and (3) the particles have negligible effect upon 
the fluid motion. 

The particle drift velocity is a maximum for certain opti- 
mum sized particles. The optimum particle radius is given by 



in the combustion chamber. 

The maximum particle drift velocity is independent of the 
chamber size although the optimum particle radius is a function 
of the acoustic frequency and hence of the chamber dimension. 

The magnitude of the particle drift velocity In a typical 
liquid propellant rocket chamber varies from . 4 to 4 feet per sec- 
ond and in an average particle or droplet life time the actual 
movement will amount to just a few tenths of an inch. 

It is concluded that, for rocket chambers of reasonable 
size, the drift movement of particles has a small, if not negli- 
gible, effect on the distribution of fuel or oxidizer in rocket com- 




v/here P is the frequency of acoustic oscillations 



bustion chambers. 
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TABLE I 

PARTICLE DRIFT VELOCITIES IN THE CROSS SECTION OF A 
CYLINDRICAL ENCLOSURE 






(degrees) 



ykr 

(ft/ sec) 



w 

(ft/ aec ) 




(ft/ sec ) 






. 1 


0 


11.8 


0 


11.8 


. 1 


22.5 


10.6 


-.3 


10.6 


. 1 


45 


7.5 


-.4 


7.6 


.1 


67.5 


4.5 


- • 2 


4.5 


. 1 


90 


3.9 


0 


3.9 


.5 


0 


37.4 


0 


37.4 


.5 


22.5 


34.4 


-5. 3 


34. 8 


.5 


45 


27. 1 


-3.7 


27.4 


.5 


67.5 


19.7 


-2.6 


19.9 


.5 


90 


16.7 


0 


16.7 


.9 


0 


10.9 


0 


10.9 


.9 


22.5 


12. 2 


-6.9 


14.0 


.9 


45 


15.4 


-9.S 


18. 2 


.9 


67.5 


18.5 


-6.9 


19.8 



90 



19.9 



19.9 



9 



■ Pressure Node Line 
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FIG. I -GAS VELOCITY MAP OF GAS OSCILLATING IN THE FIRST TRANSVERSE 
MODE IN A CYLINDRICAL ENCLOSURE 
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DRIFT VELOCITY MAP FOR PARTICLES IN A CYLINDRICAL ENCLOSURE 
CORRESPONDING TO GAS MOTION OF FIGURE I 



I 
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APPENDIX A 

SOLUTIONS TO EQUATIONS 3. 16 TO 3. 21 
In this appendix solutions are obtained for equations 3. 16 
to 3,21 which are developed in Chapter III, 

If a particle is considered at rest at CT= o at a point Z , , 

then the initial conditions for the equations are: 



= Z ^ (p^(o) = (p ^ = 0 

2,(0) - (0) = (p, (o) - (o) ~ 0 

= (pjo) - ^(0) = O 



( 1 ) 






Zero Order Equations. The zero order equations, 3. 16 and 3, 19 
are: 






0 



( 2 ) 



(3) 



= constant = 2 and = constant = satisfy equations 2 and 
3 and will be used as the solution. 

First Order Equations. The first order equations reduce to the fol 
lowing upon substituting Z and into 3« 17 and 3o20: 

^ y component: 

“ a'J 






j; (ircxi) sin ? Sin 



(4) 
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L ^ component: 

. /. J ?, _ 



trey’ 






. ir Pa X, 



(5) 



The solution for (p^ is obtained first and then the solution for 
is found. 



If: 



D = 



S'* 



S/'n (p 



the component equation may be rewritten as 







D sin 



^Tr\)a. T 
C 



( 6 ) 



(7) 



The solution to 7 is determined by the usual technique of finding the 
complementary and a particular solution and adding the two. The 
constants appearing in these solutions are then eliminated by use 
of the initial conditions. 

The complementary solution to equation 7 is: 






C, + C; 



- \<y 



For the particular solution try. 



4> . 



'Try> a y n ^ ^ 
C? Sm — -f- CoS 



<5?rpa J 



(8) 

(9) 



In this appendix, the subscript C will pertain to the comple- 



mentary solution of the equation and the subscript p will pertain to 
the particular solution. 




The first and second derivatives of equation 9 are; 

pTc>a/ n ^ .arpaT j r JTpa'y 

G3 cos u -~c 



( 10 ) 



-- 



- (^-^J 






si., cos 



ay pg 3 " 
'a 



( 11 ) 
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Substituting 10 and 11 into equation 7 results in; 



/a7rva\ ^ . 

" \ C ) ^3 


a TT p 0. Q-' /37rpa\^p a.TT^ay' 

C 1 C // ^4- c 




^ 1/ s ini a- . i/a rpo. p 5 ,^ aTrpa.3' - 
c C ^ C W C 


Q sin -STTpaO" 

0 


f 

. ( 12 ) 


Equating coefficients of like terms gives: 






- ( C 3 


Q = D , 




(13) 


and 








K C, 


-i^)% = 0 . 




(14) 



The solution of equations 13 and 14 is: 



C - - ° 






- I< D 






(15) 



The complete solution for c{) is then: 



4>. = 



C, - c 



- /cT 



. a?r^’4.T 

~D Sin — ^ 

P * ■ 



KD cos 



^TTpa T 






(16) 



To determine C, and the initial conditions are usedo 



r 



(p, (0) = O = 



l<D 






(17) 



(p/(o) = 0 = 



-KC^ - 



D 



(18) 



The solution of 17 and 18 is: 
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(19) 



c. = 



D 



aryg 



KD 



I , ^"[/<^ (^yj 



The complete solution for is then: 



<>,= 



£rpg 

C 



K D 



W 



jTt^g 



_ kjr 



^ . ^Try’CL'J 

U sin 






_ X D cos 



^TrPa.y 

o 



l<^ + 



(20) 



A similar process to that carried out above leads to the solution of 
from equation 5. 



If E = 



TTo' 






( fTo^?) cos - Jo ( r^2-} CQS ] 



( 21 ) 



then equation 5 reduces to: 






f s/n 
^ c 



( 22 ) 



Since this is exactly the form of equation 7 and since the initial condi- 
tions on <f>^ and F, are identical, the solution can be written down 
immediately by analogy and it is given by: 



I 
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z. = 



£ |< E 



f- Alt)) a, ^ ~ K y 

t - 5 — e 



(_ ^ I • _ — 

(^j'J k[i<'. (i^yj 



£ sin 



jTTpg.'T 



[k% 



L' C - c 

K E cos —3 — 



<3 7rpq. 
C 



[k%(^-F"f] 



( 23 ) 



Second Order Equations. The second order component equation, 
3. 18, reduces to: 






T, cos 4? 

.. 



^ ) COS pj sin -i- (f)^ 5 in ^ 



^ T, 5/^4 ^ j -g 



(24) 



when zero order solutions are substitutedo 
If 



and 



f = 



G - 



(yo<fl,i'rr°<^)cosp - ^ Tj(7r<^?) ocs 4 



] 



r 



^ Jo('rr^?)5inp ^ 



(25) 

(26) 



then equation 24 may be expressed as: 



(27) 



+ k ^ = i:, F S/'n i2p~ tip. 6 

The values for , 7 , and (p can be obtained from the solu- 

tions given in equations 20 and 23. Thus, from equation 20: 
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~l 



J-J 



D 



D Co5 ^ 4- k O sin 

' O / (L C ^ 









idyj 



^ sTfy^:^ Kr 



c -^o^mpy 



[k’^ii^yj 






(28) 



- 3K o' S'" TT^ C »s t /<"D" S/W= 



n 






(29) 



Substituting 29, 20 and 23 into equation 27 gives; 



Y- ly _ 



r 






a ^. _ X k r_ ^ ^ (l^ye ' ^Cc b »L^^3r J' 



j z d" e~'^\irx. ^ ~y F D^cr i"" A (g^±/ 



- rZFK 0^ COS -f. B K^O^ S/^^ 

Cy c a, Cy 



[k^* pp-^yj" 



+ E F sir, lipj ' ^ ki±jm}pp2. 
>. [ k = + (i4f> y] k (pp/J 



k[K’■^■ (ZOi3-/j 



- 06P^‘^e-''y,n^P^^ + 0S4?*-‘'i«i^“^ 






k[k‘<-(ipp] 



I 



i 
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• i'Try'^^'y ,, rs/ 

k£Fu 6 S — f-- b/n -3 K D(r Sin 



0 6 S/W 



s vi^'x y 















_ K 06 OOS — ^ 



j?TK<:lT aff Wa J' 



s. T>^<y 
a. 



[k-^ (ijny-rj 



( 30 ) 



Equation 30 can be simplified by using the following trigonometric 



identities: 



5 in 



iSC^ = k %,rt± 2 ^ 



Sin 

o 



- / b _ 0 - , ^'frpa T \ 

- ^ 0^5 —3 j 



(31) 



Cos 



* - 






^ ^ ^ c« 



) 



The results of substituting the above identities into equation 30 is: 







^£Tp 5 j 


I'e' 







jy^ jy 


L 1 


[k'* 


1^ 



^ (S^;V H- 



-, -^rpa-T 2/ 2/ ,4()rua r 

H KD { ~ JS/n — +<£^0 — K D CO s ^2—3 — 

r I . ;j / D -TTlJCL 1 



^ (F£ 4CS)(^^)5'//?i^-‘"^ ^ (Ffc -hOe) k sin 

[/cv 

- (F£ s/^ ^ (f-£ -^06) ^ (f£v-PS)6ci 

k[ k UT^^l^WyT~ 
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(Fl +Dc) 









( 32 ) 



Equation 32 is the linear second order differential equation which 
must be solved. The solution is obtained in the same manner as 
was used in the previous section to obtain values of the first order 
functions and Z, • 

The solution of the complementary function is: 



2 .-C =[C, + C. 



(33) 



and the trial particular solution is: 






C,y+C.e +c^e +C^e sini2^iS 



+ C, + Cs S /a tlHaJ" ^ c, ^ C,,s,n i 



~) 



fryayy 



(34) 



The first and second derivatives of are: 












C- 



^4rya r yir)}a. r </r> + < 5 J 213 - r. c,o.^ yrya'X 



( 35 ) 






^ k % e- ). ^4- 



_x r 



■i- k 

r* 
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- Kr . 
5m 



2 Tt>>C . 

C. 

Ky 



^ - SJryL^ K C,c: ^ ^cos 



_rr 



2nyay 



'(,^ (i, " o c. 



^ / 2 ^o^ <? iT^a ^ 5//? ^ fry^ ■J ' 



4- rv’u T 
C5 



(36) 



Substituting 35 and 36 into equation 32, and equating coefficients of 
like term results in the following evaluation of constants; 



r ^ 



kc, = 



^ D' - -t (.Fi +DS) 

[■F77I|ia)7 



sk y - 






- i^S c, = 



o3 TTyfO. 



Yk‘Hs^'fr" 



k[k^* {i’^rj 



- l<(-^^)C, = 



_ / S>ry’<^\^ 






kc, = 



% + iilftrDS) 






- 122:!£. -/ii?2^V(3 :r 

C— 1 ' C^ y ^ 



^D(n) 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 






O 



(43) 



I 
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/> - 

\ C- J ^ic 



(Ifc +C4(X^ ^ 




(.-=’2''p7 





(44) 



Solving equations 37 to 44 for the constants gives: 






c, 






C4 - 






-(f£-^D&) 


1 

r r- 


>6 .os;i:^r") 1 


[i,^, 




_ 



(45) 

(46) 



c,= 






"i" (F5 D 6) 



;?*:rya 









- + \kifB+P&) 






(47) 



(48) 



= 



- f F£^ 



(i|tL?)YK%(^")7 



C,. -- 



- (F£ -hP G) 



(49) 



The complete solution for can now be written as: 



2. - 






'J 



K3' 

-t ^ 

[k u i^^rj Kih-^i- 'I ‘ 



r-,-5 /r- r ^r'-'.'/ 44.- /,- . - ^ i - A J~ SJjL?- ^ 

^ l^E Cr - (F£ -tPirjjC -j- (l- £:+!>))( c" VC — 



^ ^ ^ ^ [ii > i J [ " '"( 
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f 1k -,( ir>¥rl>-'( '¥j 1 Ir^i^ L i' ' ' ‘1 



(yt tPcr) - _ (1-t ^DG) cr:> ^ 



( 50 ) 



To determine C, and the initial conditions are substituted into 50: 



- O 



c, f + ^<^"(■‘''2'")^ ^ Ei£i_zA'-JJ±,^2L 



[- -t ;i 2 yo^ {fc ^-lg) 

L/C CC‘Ai^yj*^iK‘G*sil 



-f ^<2 { y tr t ^jg) _ 



(Ft *- ucy 






(51) 



(0 ) - D 



-k c- 

k [y^ y k [ F y (k^yj " 



, k[;pZD^-(F..P&)J ^ 

'[KU(iipy]^ [kkyj^yyfkkc^)^ 



^ k (tt +Dio ) _ ( F £ -f U G ) 

y t ( k^yj[y7(i^^ kfy^^ 

t 

Solving for C and C > results in: 

•2 ^ 



( 52 ) 



c. = 



y-^yy y^Kk yj^y] I "]■" 
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-h 



(Ff+DS) 1-^-1“^+ fK ^ --% tc VJ +^(re.tDSl _ , 



_ ^Kh-H DG ) 

K=[K=;(i!prj_ 



( 53 ) 



c,= 



-_[4 -■iim>-D6)] -i ~£D^ -|^‘'yV ((ypfc'I 

K='[k^+ °+(a|wi^j'-‘ [ 






+ 



('F£ 4 Dfe) _ § 0 ^ _ [?gp^- (Ft +£'&)] ^ ikt^D^X 

k\k\^y] 2KYk’4(tt%/ 



- L-§ (--F)V 4-ai k’D^ip] 4 (? lE'?’) 

[ k= 4 (J»)^]\4f-*)[K"4 (t£V4«J ik^-.(ia=^y]iT[kM«^^y] 



- *& (Fc 4 Co ) 



(54) 



The complete solution for 2^ is obtained by substitution of equations 
53 and 54 into 50 and is; 




- [%D“-^(FE4CS)]_ ^ [JZ0^- ('FF<T£J 



- (Fg-tCG) (igSr) =' 



_[ ? ( -HVi' - .H (l£-* ? ” I 






4 
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(Vt + OQ>) ? D P^- + 0(-'^)J 



- r- % (■^*'-^) tj?K=D" J 

_ ( Ffc + £^'L'') -f (('£■+ L 'S') -*■ i- j,^ 






K J- 



72iy^ \ ^ 



. fi i D' - ( ht + " ) JC “ -j. (Ptr -?-C/'j)('^'^) 






. ? -i 



./cr 



*<■ t-' ~ ^ ^ A (jf 



-f-|^2 X- -i l/'tr + CL'^J X 

Ki 



'2 0^ ( ^ ^ 



-f [ ^ ^ 0"'" (^tv-Okj ^6 

[i<Mii‘>/] 



-t.3.,. ;r/« J- (tttt,..y(-<5'';t'\’',^ "-^r---'' 






_J 



_, j_- ?6 r-j"' )'D^‘ + •<^>:-y±Tl 



^'Hk‘f (if-'V’J ^T'l k v(lf>)'J [ k " n'-r )'7"( i;"“J( » C-' "■•' ' ) L; 



^ Ifyii 






I, . (ff fL.n/^'-d 



y 






( 55 ) 





I 




« 
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The second order component equation can be treated 

in a similar manner. Thus, substitution of zero order solutions in 
3. 21 leads to: 



dj JJ Jo Jo Jj 



If 



f~l = sin ^ Jt> ( )_J 



(57) 



and 



I = [n;(!r-'2)cci.f ] 



(58) 



then equation 56 can be expressed as: 






c/J 



Tj 



[yi, H b/h f-c^^ J 5y'ri 



(59) 



Terms involving Z, and 4>, are evaluated by use of equations 20 and 
23 and are given by: 



Jj, _ 

jy 



- 



' o 



I CsS + K t 5/n 

'■ 



(60) 



US^JWY] 



(61) 
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J7 J7 



—Z' 



21 S. 



[K^^ (i^yj 



+ 2D(i^)KE g-^^sin -f- (£|V£^J^D£ 



m^v] 



21 ^ 



sm co^ -f- k^DB ^/n^ 






[ic-^ {2j^yi 



(62) 






- K'D (^/‘^)‘^D5in^'''^-, 

- c — 1 ( 63 ) 


















+ £ (-J^TO Sin 
k[K^ 








[k^f(i^^Y]^'~ 



-(i?Z<?) /C£D b/n ‘i?:^^coi iJTk’^J ^ k^tD fw«) g- 



[k^(E«»)‘j^ 









(~) ktc, 5« - (2-p')x"£L i*7-"- '>' 



2Vy’CL 



Lk^-^i^-^-^y]"- 



(64) 
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2 = 
'j:r ^ 



^ Q + KED^^- sin 

--- ki?tJ^WJ^ 






-EC Y^-^'' V Y-^Ye oe '"'ivs - K ED e' Yin srE^r 



— KJ' ^7»'‘^^» ^ 
P 1 «i- 



,XT 



[EY^Y] 



-£DY-Y-^]e'^'^sln^J-.Y? +i^~'^^CtCC5 sla 

\ ^ __ ^ c _ a 



-KtO sin iy£L? -Kao (i->7!^) i - '■ - i ^i- 



a?-i 









■aVr^^Ktl, - KV-Ci<-^ 



--r [*■' 



c 



c 



Substituting 62, 64 and 65 into equation 59 gives: 



^1 = 
o/O" 



(65) 



2~[ k^ 



-4D(£2^)Kfe-'"^//?^^^ -^(^Yd£ 

2 pcapry]^ 



-^4(^)DKf sm - JkYe siY 

zCY^ 
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•j-S-KBD 1 


[^] 




l"Ds<^^Y“^--*f’"‘‘AiO 


zK fK-'+ 


1 ^ 






f k Lli-', 

'ii t't (i^r'-';"^ ‘ 

+2(S.Y'=-)^^0 -M -’!/'* -’■' - ^k'‘tOi^-’:'-^IC-'i>siq ^ 

^ 

+ = (^^rkSD t 3 (i^)k’E 0 C0^‘ilp2- 

? (=^j 

-SKeD(3k^]=s- *"■+ J k (^j’fOWS 

# K fK-' ■/■ ^i^vr 

-Jk%D(^^^)£''^^-^S (^^) K^BD -SK^^O iin 

5 (i^) L 

+5/t£D (^§i!■/e-^''' j j Ve D e - ''■^c« 

? A, U-+ (izy-f]- 



j 
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^ I '■ i- (» 

-f (^S^) k £ D «s £ ^J i- SKi o S/n ^ 

2ik'r 









kl\^f (^-Yj 



+ (EH+0l)kV^if''“-'' - (t:Hnil^E'’'^t^-^\}r, 

_ Xl^ 

klKU (J ^/'^yj 



-(EH^CDiin’^S^ - ieH*Cl)Kccy-^E-' -"‘ --- c--"" I 

[k't ?-y“l. K=+(2K->)i/ j. 



Using trigonometric identities, equation 66 is simplified to the 
following: 



^ 1^ d 4.2 — 
dr"- 1 3 - 






U'V yjry ' ff 



3.{iY^)k£D c-'^^sm - il^fyy^Dk C.c: 



^ B pa y 



i- [i-y’"y] 
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fJ(£^)kBD 5 m 

Pif/r^-. z kck^^ L-^yj^ 

-^m^*)KED Sin -^^ {^-^YbO G.-^'^sin 

p [ (£J^y/ ^ KiK^t 



-:ik^BDslni^^ + (EH ^VJ)i^)s)n^:!^ 






A L i^^y] 



(£//+DJ) A i,n - (BHH)J){^^)e'^^s/r2 



J'P-PCL 

c 






A A A"V (^'p'-yj 



-h<.th-tV2)-r hi Lh^j J- 









(67) 



2 J 



The complementary solution is again: 

'f ic jpi + ^ ] 

The trial particular solution is; 



( 68 ) 



~ |Cj CT + C 4 -t- Cf e ^ cefs "+ C* c 

The derivatives of d) are: 



- i'fr\>a.y 

^ in. r 



C„ sw aia-?" 



]. 



( 69 ) 




KC.^c^'^^-KC^ e-^^c ^5 



-iysc, C-^ln 
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-f^ C, cos C^sm Q, Cos 



(70) 



c - 



4A:t rk^CrS'^'^coi +^fHCs&^"li'^^'^‘^^ 



c 



-fK^^ {if’^%t' '^cn:. 

-^ S-p kc,e- ^ if 






\ c y 



3 . x47TpaJ 

In ->— 



- C^“)'c„ S!n j 



(71) 



Substituting the above results into equation 67 and equating coeffi- 
cients of like terms we get: 



akx^ 



- kj(fcH + DI ) ] 

[V^I- {-- J 






(72) 



(73) 






-zxtof^ - imuh^a) 

zlKH^-y]' -niK’^c-sp-yr k[K->('-r)-J 



( 74 ) 
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4DE ^ 



-a (^^)^Dl + ME Z-^ + DJ) 



( 75 ) 



[jk V j "" 



(76) 






?[k- t (- E'z- {- ? 



K'lD 



"Jl 






-{^Tyc, ^kiijyic,. = 






(77) 



(78) 



-k(^r-1C,-(q^-*)%, = 



, ; rr4<i-\-^, 



« /r>^-eV‘ 


f- 


~ c / 


J 



-f- 



(..*ci)F,r; 



^ ^ k ~^t L' -hKitzH-fL 2 J 



Solving equations IZ to 79 gives: 



Tui (79) 



C3 - 



~4 (BHi-01 )_ 

k[k^, 



-£;£: r-S" . 



(80) 



C. = 



-J0B._{-¥^ )__ - (EH i-bl ) - 



(81) 






“i pt .p'^ 



FTlfi 






4 



:, '■-i' 
^'1 


















I1 






















|i™% . j^ Pi.||^ffii ■' 
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r = -a(3 |^j"DE ^ Kt o 1 

^ 



(82) 



Cr 



a £'D _ (e-H-ypi) - 4i K‘s D 






+ 



!^| (EH+DI) 



- K'^ED 



1 






(83) 



C>-- 






+ (thi-Dj) -3nn^to 






- K'^D . ._ . ->- 



(84) 



I I r , 



^ 's k [I, ^ r-"")T - (^-I-^)’'!;'* (■=?'’“/] 






E.1^ Hj^Dl) 



-_k‘(£H + Dl) 



(85) 



C.o^ 





-E 4(^c^‘^)k£D 


■1- SK^tC 


_2Hk“+(-=r-^l7" 


1 





. (chK)3)rT’"' 



c / 



- (t N + D I ) k 

/< »va ^ , I a /; -•< 

{- c ; i k ^ ; J 
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The complete solution for is obtained by substituting the above 

constant into equation 69 and adding the complementary solution. 



I 



-65 - 



- 






k[k’f(V^y] iK^[K‘*<S^')‘] 



-SDE (^*) C - "’cs - a kVo e- d^^-fEwoOe^JiSgS?' 












-ar-p; -E (EH.Dl){^)c- 






k [k^t 



" # (1^7 

(EHi- pj') coi^=^^ - ^ k^£ Pco:si^^^ {£H+Ol)cos '^ 

(^^)’j[ k"^(±^yj 

-k%D cos - ;) |<£0 (£|y-^}^;/, ligiy t ^(Eti^Di)sM^-^-? ' 



•f '56 (anwi) S-//7 

2 [k i (-4^";:k ■[ k^.(t^->jj 



^/^n'ya\ - 4-7rV q 4'7TPa J 

- 3 ("=^ 2 r J/^tD Sm^-ZT' - k ED bin - 






?K LK"+ j ^ 2 [k^-^ 2 



- f£ ^ Di. c.>. h>DI, cr. i -£a s,r^ ^7-“’-^ 






^ 3 

J J 



-h ^U-r'")kLD ,ir^ -f- D :./n 



( 




t 
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(EH ^DI)K 



( 87 ) 



To determine and the initial conditions are used: 



<t‘J0) = 0 = 



r + - - -.-DcJzPjl-idElL±£lL 






_~h(£H I) -_ ±’'j= c _ 



^[£'^(^-¥-)Th'^(^-T-y] 



f 2 EL 



i- (Eh^PIj 7 



-i- Jk^BD + 3 k'^bl 



- 1<^(£H-^DI) 



i{£T)%'^¥¥Pf (-P'TIk^H-^-FPT 



( 88 ) 



(p[{0) = C = 



j(^ Q - ^ H ■*• plj ~*' 4 /C'* fc 

Bk[p.(£p’^y¥ £[K^.(£2r.>^yf 






_ -J.^)k\tD__ 



,_t>? ^(^r"y kLD__ 

^ ^ ("-r'’ )'j ^ 



( 89 ) 



The constants Cj and C, found from 88 and 89 and are; 



- Jji(BH +Dl) + 3Pg H- 4 k"g P 

k^[k^^(^pL^y] fkZk^^(2^yy 
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Z[k-.(i^yf[l<^(tXpr] [K%(i:p^yp^(tJ^yj 



- k^BD 



-t (EH-tDl) 



■^^{^^VbD 1 






( 90 ) 



c = 



^ (' - 4i^^EP 

k’[K^t(i^-^y] 



+ ZlC{^^)^ - J^(i:h^DI) 






*=^ 

ri- 






-j. (±^>’‘^)K^tD 



— ( h h + C 1 } 









+ DE + ^0£ + 2(EH^DJ) - aji^YBD 

2[k^^(3^yf [ic^.(^yj^ i[k^^(im)fi^\(iiks)^ 



_+_'^ 

[kk(i^y][kk(±^‘^r] 



-t^a K^EO 
^[k^^(iJ^)f[k %(if^)j 



- (^£/^+Dj) 



+ D _ 

v(%’ v[fc^- 



~ p 



— £? A ‘^iE £ 

The solution for 







( 91 ) 



is then obtained from 87, 90 and 91 and is: 



-68 = 



~ ?[;^V(i^^/j"|;i< V(irpa^j 

J (^PV^) £D -_.kL^hj±ii. 



-( fH + Ll ; ^‘^- yry'*)yD 
[k ?k Xk V ^ik y ^ 



-t2(EH ^DJ) - o2 (^T ED 



-^ JyUH + DX ) 



[k\ ?rk V [k*. (£^yj[k^(^)3 






- K ^(£ ktDJ ) 



+ D 



z [ k^f (imyj "[ kM ^ f ] /c 



~J£Di- ^-J -2K^ pC __ -JN ED -i- k Ub HfL'l ) 



-- i ( ciM - cj t- JCc e ‘*-_^' -<- 4 -r £ L c 









-^ED - 3 yj ^^ W°)£:D C ~ 

?[K ^ ( i ^^ 3 ' 7 [ K % ( ip )'] fK % V (± P ‘)'] ^ " ] 



-(•£ 52 ^ K-ec e "^' 



+ (Eh +-EJ / r +(? l £ X "£ 0 ^ 



a > - K -" 






- >^(E h ^ ODd ' 



'w 



K : h - 



- DE ^ r U 2 TJJ. 



zk \ k ' y^^n 



- JiLti^jC cos 

^ LK^^ my-fV 
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^-)3 






3 1 ^ 






S/^ c^ 



niZ 






is 



- U kB oyj^} t ~ ^ 1^ f- ^ ^fc L "i 



i” £<,'V 






(BJ i- D I) u>5 -hk'kDu^ l£'^_ 



- K ^FO Ci’s ±1^2^ 









+ ^ (Ek+Dl) sin 12p^ 

(^)[K^ (i?t!j5f*S(t¥-“/J 



- A?^f O sm 



-3(^^Vk£D S/n 

2[k’^(^-)T [k'- i^F-yj 

y F;i Lb H -i- Dt) S/n 



i^-^)[k^^i=BBB^)f[k^H±^-f] i^-T) [k^H-F-'fj 



-hJKE Dm^^fccs 4 k% 1 5 i ^yk 


-h ^ k^t L —‘■'~i^ ^yy 






^ J-YvY^ (^-/j 


-{EH-^DI) Co^ 


— /C {th^L^ / ^C :> - 




■ [H=t(kjyHyy 


(-%-//* (4”/j " 





_l_ 4yi|V-» j l^cDSini f- ^k^to 3/n ^ -{B Hi- L DL^I^’V 5>n i2LB^ S 
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_ (EHtDI)k 



The terms proportional to time in equations 55 cind 92 are 



circled and are the terms used to find the particle drift velocity o 




# 
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APPENDIX B 

A DETERMINATION OF THE JvlAGNITUDE OF THE FIRST ORDER 
DECAYING EXPONENTL^L TERMS OF THE PARTICLE VELOCITY 
The first order decaying exponential terms iii the solution 
for the particle position in the cylindrical enclosure are: 

_ lt"T' 

c j€ 



2,= 






K[K' 






TJ 






Here, 



0 ~ sin $ 



( 1 ) 



£- = ir^ cosp - 3;(7rc<,^z)] 

k- =r ^ 'Tr'J^ <=T' CL 
r7n c, 

From equation 4. 1 the expression for the particle velocity, contain- 
ing just the first order decaying terms, can be written as: 



r\r = 






CS 2^- CS £ 

dr ^ syj 



S dA 1 



( 2 ) 



L 



Substituting the expressions for ^ and into equation 2 the fol- 
lowing is obtained 



nr = 












bed D\-q ; c 






(e\ 



u / 



(3) 



For an initial position of the particle of £ = ^,1- | = 45®, for 









= . 316 2 uid for gas and particle properties evaluated for the 



I 



! 



-72- 



conditions used in the calculations for Table I, the following ex- 
pression is obtained for the particle velocity due to the damped 
terms: 



_ / .^35oT\ 

nr= 4. (7/./ ; 
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_ ^350t\ 

^ j 



( 4 ) 



For "t = . 0023 seconds, which is time for approximately two os- 

-^Z5oT _,q 

dilations in the chamber, (r — e , and the velocity amounts 
to: 



or 



C (-• ooz) -h 2} C Oo4) 



fV 

'"^sec. 



( 5 ) 



For these same conditions the steady state drift velocity as calcu- 
lated in Table I is: 

= t (Z7) (-.^1) ^y^ec. ( 6 ) 



Comparing equations 5 and 6, it is seen that after just two oscil- 
lations the velocity contribution of the decaying term is negligible 
compared to the drift velocity. 



I 




r 






I 








I 



I 

i 



